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We compute the genuine SUSY one-loop quantum contributions to flavour-changing MSSM 
Higgs-boson decays into 6s and sb using the full diagrammatic approach that is valid for all 
tan f5 values and do not rely on the mass-insertion approximation for the characteristic flavour- 
changing parameter. We analyze in full detail the dependence of these flavour-changing partial 
widths on all the relevant MSSM parameters and also study the non-decoupling behaviour of 
these widths with the SUSY mass parameters. We find that these contributions are sizable as 
compared to the SM ones, and can be very efficient as an indirect method in the future search 
for Supersymmetry. 



1 Introduction 



Rare processes involving Flavour Changing Neutral Currents (FCNC), and in particular, the ones 
related to Higgs physics, provide an extremely useful tool to investigate new physics beyond the 
Standard Model (SM). In the SM, there is a strong suppression of these FCNC due to the GIM- 
cancellation mechanism. In contrast, the Minimal Supersymmetric Standard Model (MSSM) 
provides a natural framework where induced FC scalar interactions could be significant if the 
soft SUSY-breaking mass terms contain some non-diagonal structure in flavour space. 

Here, we are going to concentrate on the neutral MSSM Higgs-boson decays into bs and sb 
that are induced by SUSY loop contributions, in a general scenario that includes both possible 
flavour changing sources, via the Cabibbo-Kobayashi-Maskawa (CKM) matrix and via the mis- 
alignment between the quark and squark sectors. We find that these FCNC effects are large, 
considerably enhanced with respect to the SM contribution, and therefore, they are sensitive to 
search for indirect SUSY signals. We also checked that our results are perfectly in agreement 
with the ones of the effective-Lagrangian approach in the large tan (5 and small FC parameter 
A limit. 



One of the most important features of these observables are the remaining non-vanishing 
contributions even in the limit of very heavy SUSY particles and, in addition, that they are 
enhanced by large tan/3 factors. Such a non-decoupling behaviour of SUSY particles in the 
Flavour Changing Higgs Decays (FCHD) can be of special interest for indirect SUSY searches 
at future colliders, as the forthcoming LHC and a next e + e~ linear collider, in particular if the 
SUSY particles turn out to be too heavy to be produced directly. 

Here we performed an exact computation of the complete SUSY-QCD and SUSY-EW one- 
loop contributions from squark-gluino, squark-chargino and squark-neutralino loops to the 
flavour-nondiagonal decay rates of the three neutral MSSM Higgs bosons, and therefore, our 
computation is valid for all values of the characteristic parameter measuring the squark-mixing 
strength A and for all tan (3 values. 

2 Flavour-changing interactions in the MSSM 

In the MSSM there are two sources of FC phenomena. The first one is common to the Standard 
Model case and is due to mixing in the quark sector. It is produced by the different rotation in the 
d- and M-quark sectors, and its strength is driven by the off-diagonal CKM-matrix elements. This 
mixing produces FC electroweak interaction terms involving charged currents and, in particular, 
in this case, SUSY FC electroweak interaction terms of the chargino-quark-squark type. The 
second source of FC phenomena is due to the possible misalignment between the rotations that 
diagonalize the quark and squark sectors. When the squark-mass matrix is expressed in the 
basis where the squark fields are parallel to the physical quarks (the super-CKM basis), it is in 
general non-diagonal in flavour space. This quark-squark misalignment produces new FC terms 
in neutral-current as well as in charged-current interactions. 

Assuming that the non-CKM squark mixing is significant only for transitions between the 
third- and second-generation squarks, and that there is only LL mixing, the squark squared-mass 
matrix that has to be diagonalized via a 4 x 4 matrix, R^ d \ in the (sl,s_r,6l,6k) basis (similarly 
for the (cL,CR,tL,ti{) basis with rotation matrix R^), can be written as follows, 

M 2 Ls m s X s \ LL M L , s M Lfi 

m s X s M 2 Rs 

\ LL M L>s M L>h M 2 b m b X b 1 ' 

m b X b M Rb t 

where M 2 , M R and m q X q are the usual flavour preserving entries of the MSSM squark squared 
mass matrices. The previous assumption is supported on many mSUGRA inspired models. 

In our parametrization of flavour mixing in the squark sector, there is only one free pa- 
rameter, \ll = A, that characterizes the flavour-mixing strength. Obviously, the choice A = 
represents the case of zero flavour mixing via misalignment. 

After performing the whole diagonalization, we computed the loop-induced flavour-changing 
neutral Higgs boson decays into second and third generation quarks whose detailed results can 
be found in 1*^=1 In the next section we will discuss the dependence of the decay rates and 
branching ratios on the MSSM parameters and A. 

3 Numerical analysis 

Here we numerically estimate the size of the loop-induced FCHD as a function of the MSSM 
parameters and the mixing parameter A. The GUT relations M3 = a s /a Sy^M^ and M\ = 
5/3 s 2 ^ /cyy M2 are assumed. For the numerical analysis of the FCHD rates, only values of A (in 
the range < A < 1) that lead to physical squark masses above 150 GeV will be considered. 
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Figure 1: SUSY-QCD F{h -» bs + sb) (GeV) as a function of (m A (GeV), tan/3), (M (GeV), A (GeV)) and {fx 
(GeV), tan/3). Fixed values chosen are n = 1500 GeV, M = 600 GeV, Mg = 300 GeV, A=200 GeV, m A = 250 

GeV, tan (3 = 35 and A = 0.5. 
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Figure 2: SUSY-EW F(ho — » bs + s6) in GeV as a function of /x (GeV). Triangles, boxes, and diamonds correspond 
to m A = 200, 250, 400 GeV, respectively (left).r(/i -» bs + sb) as a function of M (right). 

The MSSM parameters needed to determine the partial widths T(H X — ► bs + sb), for H x = 
h°,H°,A°, are the following six, tua, tan/3, //, Mi, Mq, and A, where we have chosen, for 
simplicity, Mq as a common value for the soft SUSY-breaking squark mass parameters, and all 
the various trilinear parameters, A, to be universal. These parameters were varied over a broad 
range, and grouped into different pairs in order to visualize the individual dependences of the 
FCHD widths for each neutral Higgs boson. The most important thing to notice is that we 
obtain large values of the decay widths. Besides, as expected, the SUSY-QCD contributions 
dominate by large the SUSY-EW ones, being at least one order of magnitude larger. On the 
other hand, a common clear behaviour of all three decay widths is the increase with tan/3, 
yielding maximal FC effects at large tan/3 values. In the following we comment shortly on the 
ho and Hq decay rates. For the Aq, similar results (not included here) than those for the Hq, 
are found. 

In fig. E] we show the behaviour of the SUSY-QCD contributions to T(Jiq — > bs + sb) with 
the different MSSM parameters. We can appreciate the above mentioned grow with tan (3, the 
expected decrease with Mq, an approximate symmetric behaviour under \x — ► —fx, a growing 
with \fj,\ for moderate \i < 600 GeV, and a nearly independence with the trilinear parameter A. 

We show in Fig. [21 (left) the behaviour of the SUSY-EW contributions to T(ho — ► bs + sb) 
as a function of the fi parameter for three different values of tua- The shaded regions in these 
figures correspond to the region excluded by LEP bounds on the chargino mass \fj,\ < 90 GeV. 
The width for the H° decay (not shown) is approximately symmetric under \i — > —\x, depending 
of the iriA values, while, as we can see, the T(h° -^bs + sb) width is more asymmetric with 
respect to the sign of [x, but all decay widths increase with for < 500 GeV, then reach 
a maximum, and finally decrease. Regarding the behaviour at very small /x values, we have 
also found that the widths do not vanish at \i = 0. The origin of this comes entirely from 
contributions driven by electroweak gauge couplings. 

Fig. |21 (right) shows the behaviour of the SUSY-EW h decay as a function of Mq. The h° 
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Figure 3: SUSY-EW (Diamonds), SUSY-QCD (stars) and total (boxes) contributions to the Br(H x -> bs + sb) 
H x = (h°,H°) (a) and (b). The prediction for Br(6 — * sj) in the MSSM is shown in (c) as a function of A for the 
chosen MSSM parameters specified in the text (solid line). Short-dashed horizontal lines limit the experimental 
allowed region and long-dashed ones represent the spread out in the MSSM prediction induced by the theoretical 

uncertainty in the SM contribution of 99 % C.L . 



decay width has a small value for light M$ due to the fact that chargino and neutralino contri- 
butions have opposite sign there. For higher values of Mq, the neutralino contributions change 
sign and the partial cancellation disappears, therefore, the decay width increases until it reaches 
a maximum and then decreases for heavier squarks. The previously mentioned cancellation for 
small values of Mq is less obvious for the heavy Higgs, the clearly visible effect is the decrease 
due to the growing squark masses, which is slower in the latter case. 

In the following we study the behaviour of the corresponding FCHD with respect to A. 
Investigating the contributions from gluinos, charginos and neutralinos separately, we found 
that the gluino/neutralino contributions increase monotonically with A, being exactly zero for 
A = 0, as expected. On the other hand, the contributions from charginos show explicitly the 
two FC sources, CKM and quark-squark misalignment. In fact, the chargino contribution is 
different from zero for A = 0. The non-zero value at A = is due to CKM mixing which is not 
present in neutralino (or gluino) loops. The CKM effect and the effect of squark mixing compete 
in the SUSY-EW case in some regions of parameters but for larger values of A the non-CKM 
flavour-mixing effect dominates. 

In fig. (a) and (b) we can see the SUSY-QCD, SUSY-EW and total contributions for 
m A = 400 GeV, tan/3 = 35, M Q = 800 GeV, A = 500 GeV and M 2 = 300 GeV. Note that 
the absolute value for gluino, chargino and neutralino contributions grow with A separately and 
that the decay rates of h are much larger for smaller values of uia and therefore yield larger 
values of the branching ratio Br(hP — » bs + sb). Notice also that not all values of A plotted in 
figs. 01 are compatible with b — ► s~f constrains for the chosen values of the MSSM parameters 
(see fig. 01 (c)). However, we can still find values of A compatible with b — > sj that still produce 
large Higgs branching ratios. For instance, the value A = 0.3 is allowed, and the branching ratio 
is around 2 x 10 -4 for ho and 10 -2 for Hq, which are several orders of magnitude larger than the 
SM value, Br(HsM — ► bs + sb) ~ 4 x 10~ 8 for rriH SM = 114 GeV. A more exhaustive analysis 
on the allowed values for the SUSY-QCD corrections of these MSSM Higgs boson branching 
ratios being restricted from their correlation with the radiative B-meson decays (b — > 57), can 
be found in'"'. 



4 Non-decoupling behaviour of heavy SUSY particles 

One important feature of these observables is their non-decoupling behaviour, which means 
that the FC effects remain non-vanishing even in the most pessimistic scenario of a very heavy 
SUSY spectrum. The origin of this non-decoupling behaviour in these contributions is the fact 
that the mass suppression induced by the heavy-particle propagators is compensated by the 
mass parameter factors coming from the interaction vertices, this being generic in Higgs-boson 
physics. The non-decoupling contributions to effective FC Higgs Yukawa couplings to quarks 
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Figure 4: Non-decoupling behaviour of SUSY-EW T(h -> 6s + s6) in GeV with M = /i = A = M- g = Ms, 
Mi = %%M- g and M 2 = ^M s , and for tan/3 = 35, A = 0.5, m A = m H o = 250 GeV and m h o = 135 GeV. Exact 
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one-loop results in solid lines and expansions in dashed lines (left). SUSY-EW V(ho — » 6s + sb) as a function of 

A (right). 



have also been studied in the effective-Lagrangian approach EEI. We will use here instead the 
full diagrammatic approach which has the advantage of taking into account all SUSY loop 
contributions (not just the dominant ones) and is valid for all tan f3 values. Since, on the other 
hand, we are not using the mass-insertion approximation, our results are more general, being 
valid for all values of the FC parameter A. We also checked that our results converge in the large 
tan (3 limit and for small A values to the mass insertion approximation results of the effective 
Lagrangian approach. 

In the following we present the results for the form factors Fl,r defined as 
iF = —igu q [Fl q (H)Pu + F R q (H)Pji]v q >H in the large Msusy limit, keeping just the leading 
^( iUgfrTv )° term in this expansion and where we assume a common Msusy mass scale, 



m b 



' 6-7T 2m\v cos f3 
,(x) _ Qew rob 



(4 X) +tanpa ( i x) *)F(X) 



An 2mw cos (5 



1 



8m?*/ sin 2 f3 



V th V cs m 2 + V cb V ts ~~ 2 



rcb 

' CKM " CRM' 



mt)F(X) 



, | ycb yes 2 , y 
r I v CKM V CKM UL C ' V ' 



tb 



rV, 



ts 



CKM v CKM m t 



J(A) 



1 r 

4 



(ycb 



< KM^CKM + ^CKM^CKM ) "^(^) 



+ 



(C^m + V' b KM V? KM ) F(X)] ] (4 X) + tan (3 a[ x 



[■)* 



(2) 



F, 



(x) 



OiEW 



L-o 



Att 2mw cos j3 



1 + - tan 2 
9 



9w) (cr ( 2 x) +tan/?<r{ 



(x)* 



F(X) 



(3) 



where, F(X) = + 1) In (A + 1) + (A - 1) ln(l - A) - 2A], J(A) = ^[(A + 1) ln(A + 1) - (A- 

1) ln(l — A)], o - ! = (sin a, — cos a, i sin j3) and = (cos a, sin a, — i cos (3). The results for Fr 
are like the previous ones but replacing m b — > m s , m c — > m( and taking the complex conjugate. 

In the following, we show graphically some of the main features. For brevity, we choose to 
show in fig. |1] just the SUSY-EW contributions but the SUSY-QCD ones behave in a similar 
way. We can see that, for large values of Ms = Msusy, the exact partial widths tend to a non- 
vanishing value, characteristic of the non-decoupling behaviour. Besides, the exact results are 
very well described, even for moderate Ms values (say > 600 GeV), by our previous asymptotic 
results, which make these short formulas very useful for future phenomenological Higgs boson 
studies. 

The dependence of the total chargino and neutralino contributions on the FC parameter A 
are shown in Fig. 0] (right). The exact one-loop chargino (neutralino) contribution is plotted 
in solid (dotted) lines and the dashed lines correspond to their approximate asymptotic results, 
given by the previous formulas. The lower lines take into account both possible FC effects 



while the upper ones correspond to the charginos contribution setting Vckm = I- We can 
appreciate that these effects interfere destructively. We can also see that the neutralino and 
chargino contributions depend strongly on the chosen parameters and can be comparable in size 
for a large range in A, therefore, neglecting the neutralino contribution, as usually done in the 
literature, is not a good approximation. 

Regarding the so-called decoupling limit where tua » rriEw, we see clearly that the SUSY 
loop contributions, in the ho case, go to zero, recovering, as expected, the SM result. 



5 Conclusions 



We have computed the genuine SUSY one-loop quantum effects to flavour-changing MSSM 
Higgs-boson decays into second and third generation quarks using the full diagrammatic ap- 
proach and therefore our results are valid for all tan (5 values and for all values of the flavour- 
mixing parameter A. We find that for moderate tan f3 and A values, some of the contributions 
usually neglected in the effective-Lagrangian approach can be sizable and, therefore, a realistic 
estimate of the branching ratios should rely better on the full diagrammatic approach. After 
analyzing in full detail the dependence of the FCHD partial widths, with all the relevant MSSM 
parameters and A, we found large rates and these are very sensitive to tan/3, fi, and A. The 
branching ratios grow with both tan (3 and A and reach quite sizable values in comparison with 
the SM ones, in the large tan (3 and A region. For instance, for the following choice allowed 
by the b -> 57 constrains, A = 0.3, tan j3 = 35, m A = 400, M Q = 800 GeV, A = 500 GeV 
and M 2 = 300 GeV, we found branching ratios of O(10" 4 ) for the h° and O(10" 2 ) for the A 
and H°, which are both some orders of magnitude larger than Br$M ~ 0(1O -8 ). As expected, 
the SUSY-EW contributions are subdominant with respect to the SUSY-QCD ones, but they 
contribute with opposite sign and important interference effects, which modify the SUSY-QCD 
effects remarkably, appear. 

An interesting feature of these radiative contributions is their non-decoupling behaviour for 
large values of the SUSY particle masses, i.e., even in the most pessimistic scenario of very 
large SUSY mass parameters, FC effec ts r emains and they can be sizable (a set of analytical 
asymptotic results can be found in refs.^21). 

In conclusion, the results presented in this talk indicate that FCHD constitute an interesting 
scenario for indirect searches of supersymmetry, with important contributions in some regions 
of the MSSM parameter space, that remain even for a very heavy SUSY spectra. 
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